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STUDY  I'O  DKVKLOP  GUADlONffiTKK  TKCHNHjUI':S 


l_.  INTKUDUCTION 

I 

There  are  three  movliit;  base  (ji'avity  ^radiometers  currently  under 
devi'lopment . The  In.struments  are  being  developed  at  Hughes  Research 
Labs  [Ref.  1],  the  Hi' 1 1 Aerospace  Corp.  [2],  and  the  Charles  Stark 
Draper  Lab.  [S,  1].  Tlie  design  goal  for  each  of  the  sensors  is  I Eotvbs. 

Kurthermore , 0.1  Eiit  vb's  (E)  accuracy  should  be  feasible  from  an  orbiting 

gravity  gradiometer  [31.  Thi'  group  of  instruments  includes  sensors  ; 

designed  specifically  to  nu'asure  tlie  gravity  gradient,  as  well  as  sensors,  j 

i 

which  utilir.e  existing  accelerometers  to  provide  a gradient  estimate.  ’ 

The  Hughes  and  Hell  instruments  rotate,  thus  modulating  the  information.  j 

I'his  rotation  transfers  the  gravity  gradient  signal  to  a higher  f j 
quency,  quieti'r  part  of  tlie  spectrum,  and  can  separate  the  signal  from  i 

some  sources  of  instrument  bias.  Tlie  Draper  Lab  sensor  measures  the  i 

gradient  signal  at  /.eio  frequency  and  uses  a sophisticated  flotation  j 

suspension  system  to  isolate  the  sensing  element  from  errors  induced  | 

by  rotation  and  jitter.  A system  of  at  least  throe  instruments  of  any 
one  type  is  required  in  order  to  provide  a complete  gravity  gradient 

5 

tensor  estimate.  Often  times,  however,  it  is  possible  to  extract  all 
the  information  that  is  required  from  a single  one  of  these  instruments. 

The  primary  objectivesof  this  paper,  oriented  toward  the  use  of  a 
workable  gravity  gradiometer  as  a sensing  element  in  several  applica- 
tions, are  given  below; 

1)  To  develop  models  for  gravity  gradient  anomalies  and  gravity 
a noma  lies, 

2)  To  evaluate  several  methods  of  on-line  instrument  bias  estima- 
t i on , 

H)  To  determine  the  performance  of  a gradiometer  in  mapping  the 
earth's  gravity  field, 

4)  To  assess  inertial  navigation  systems  augmented  with  a gravity 
gradiometer. 
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The  work  uniii'r  this  couLrucl,  and  the*  work  undt'r  a soparuto  con- 
tract (Goddard  SFC  " NAS  5-219Gt))  to  evaluate  the  performance 
of  a geodesy  mission  with  orbiting  g ladiometors  have  some  ovi'rlap, 
especially  In  the  groundwork  areas  o 1 instrument  performance  and  capa- 
bility, and  gravity  models.  As  such,  we  hope  it  will  be  valuable  to 
supply  some  results  from  work  performed  under  that  contract  in  this 
Final  Report. 
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(■HAVITV  AM)  CiKAVtTY  CKADIENT  MODELS 


sysli'ms  In  fthicli  t lu-  gravity  ^;I■acl  ionu-lor  is  an  ossontial  coraponoiU  , 
somo  modol  of  t lio  t ypi'  of  inpuLs  is  rmiuired,  Sovoral  of  tlio  models 
are; 


1)  tessera  1 linimonio  models  iif.  i n^  Kaula's  rule  [7l  , 

2)  tesseral  harmonle  models  usinn  Allan's  rule  [ !S^  , 
S)  point  mass  and  line  mass  models, 


•1  ) exper  t nun  t a 1 1 y determined  second— order  random  process  models. 


Once  a given  model  is  eiiosen,  the  gravi t at  Iona  1 potential,  force, 
and  gradient  can  l>»’  ascertained,  'flu'  clioici'  of  which  model  is  to  be 
used  in  a given  implementation  di'pends  on  the  dynamic  range  of  interest 
for  a given  system  (which  is  in  turn  dependent  on  the  system's  speed 
relative  to  the  earth).  I'his  can  vary  substantially,  fi'om  fixed  base 
application,  to  ship  spi'cd , airplane  speed,  and  finally,  orbital  spei'd. 
The  various  models  also  result  in  different  degrei’S  of  complexity, 
fins,  loo,  must  be  cotis  idi'fi'd  for  atiy  mechanization  involving  gravity 
grad lomet  ers . 


■fills  chapter  contiiins  the  expressions  for  t lu'  gravity  forces  and 
gradients,  and  their  correlation  based  on  various  gravity  models. 


A.  fESSEKAl.  IIAUMONIC  MODELS 


The  most  general  expression  of  an  arbitrary  function  over  the 
surface  of  a splure  that  satisfies  the  potential  equation  is  in  terms 
of  tessi'ral  harmonics.  flie  tesseral  harmonic  expansion  of  the  earth's 
gravitational  potential  is 
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(1) 


tthi'rt'  t li‘‘  position  of  a test  point  in  terms  of  the  spln'ricnl 

I'oord  mates , radius,  latitude,  and  loncltude;  n = GM  ; H radius 

li)  0 

of  the  earth.  P^.  (x)  are  normalized  associated  lit'nendre  polynomials, 

and  C,  and  S ai'e  normalized  tessi'i-al  harmonic  coefficients,  which 
( m I 111 

i;ive  the  amount  of  I'ai'h  harmonic  present;  the  pi’rt  urbat  ion  potential 
due  to  a particular  harmonic  is  simply 
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and  in  writing  (2),  the  phase  information  is  lost. 


(2) 


Islimates  of  the  magnitudes  of  the  J s are  given  bv  Kaula  [71 

( m 

and  Allan  [81.  Kaula's  familiar  rule  of  thumb  for  J , is 
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while  Allan's  more  eomplicati’d  formula  is 
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Although,  different  in  form,  both  formulas  agree  quite  closely  over 
huge  ranges  of  the  subscript  (,  and  with  measured  values  of  the 
tessernl  harmonic  coefficients. 

With  estimates  of  the  magnitudes  of  the  .1^^^  in  (3)  and  (-1),  it 

becomes  possible  to  estimate  the  magnitude  of  the  perturbation  potential 

\,  , and  hence  of  the  force  (?)  and  gradient  (^)  perturbations,  since 
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(5b) 


It  tun  hoi’  assumo  that  t ho  ptiasos  of  t lio  liariponics  art'  raiuioni, 
it  bocomos  possible  to  obtain  c losocl— form  approximate  expressions  for  the 
variance  of  t tie  porlurbalion  foi'oos,  the  variance  of  the  perturbation 
^travlty  ^trudlent,  and  the  covarianoo  of  tlie  pert  urluU  ion  forces  and 
uradlonts  duo  to  tlio  tossoral  harmonics. 

riio  total  variance  ol  the  force  Is  computed  for  a representative 
force  component  ((j.^  via  Kaula's  rule.  It  is  piven  by 
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tthere  fj..  is  i he  radial  force  perturbation  due  to  the  (th  liarmonic. 
The  extra  { in  (h)  is  due  to  the  fact  there  are  f distinct  harmonics 
(m  1 to  i)  for  a particular 

fj.  is  obtained  as  the  radial  dirivative 
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since 
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ttluM'o  ^ is  the  j;ravltat  ional  accoloratlon  at  the  surface  of  the  eartli. 

Squaring  (7),  inserting  it  into  (G),  using  that  (sin  q^)  I I 

f m ' avg 

by  definition,  and  using  Kaula’s  rule  for  J wo  obtain 

t m. 


An  identical  approach  leads  to  closed-form  expression  for  the  itrad- 
lents.  There 
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eventually,  usln^t  Kaula's  rule, 
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Substitution  of  Allan's  rule  (-1)  into  (01  p:ives  a slii;htly  more 
cumbersome  series  to  sum.  The  results  are 


12. 2x10 


-10 


8(0. OOl  (^2 


1.28  X 10 


U \2(4.4 


(111 


T - X 


whe  r<' 


(o.oof)  . 


I 


S 1 i 1 1 1 a I 1 V 


12. 2x10 


I.UXIO--'  ^ 

3 

(1  - X) 


The  results  ol'  Kiis.  (H),  (10),  (II),  anil  (12)  are  plotted  in  Kip;.  1. 

It  shows  the  standard  deviatioiiol  the  force  perturbations  and  pradient 
perturbations  at  various  altitudes.  Kipure  1 also  shows  that  as  the  test 
point  approaches  the  earth's  surface,  the  standard  di'vi.ition  of  the 
pradlent  usinp  Kaula's  rule  blows  up,  a phenomena  that  physically  does  not 
occur.  Althoupli  Kaula's  rule  i.s  in  common  use,  it  evidently  does  not 
attenuate  the  hiph  frequency  components  rapidly  enouph  to  produce  a finite 
standard  deviation  at  the  surface.  Allan's  rule  does  provide  for  a more 
rapid  attenuation  of  the  hiph  frequency  components,  due  to  the  presence 
of  the  (0.93)'  term  in  (T),  and  does  result  in  finite  force  and  pradient 
variances . 

One  useful  feature  of  the  c losed— form  solutions  of  the  standard  de- 
viations is  the  ability  to  determine  what  amount  ol  variance  is  due  to 
different  portions  of  the  spectrum.  A particularly  simple  example  is 
obtained  from  the  force  components  in  (11).  The  analytical  expression 
for  Oj  title  to  harmonics  above  some  depree  k is  easily  found  to  be 
1.28  X 10  ^'xl''2/i  _ X where  again  x ^ 0.93(l{^r)^.  Near  the  sur- 

face of  the  earth,  then,  the  percent  of  variance  contribution  due  to 
harmonic  degree  greater  than  or  equal  to  k Is  given  as 
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g COLUMN 
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Solid  line  is  Eotvos  units  column  or  ordinate. 

Hash  line  is  g column  or  ordinate. 

FIG.  1 1 r VALUES  FOH  FORCE  AND  GRADIENT  PERTURRAT IONS . 
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wluM-»'  X ■=»>  (().9;t)*!  In  oi'ilfr  lo  account  for  only  50'/!'  of  l lie  noise,  a 
moilcl  of  fifth  (li‘t;rcc  (~  hO  coefficients;)  Is  required,  for  75'"<,  of 
the  noise  lOtli  (-'ll!)  c oi- f f i c i en  t s ) degree,  etc.  flii'  gradient  components 
attenuate  even  less  lapidly  requiring  an  evi'n  higlu-r  order-  model.  This 
ci-rtainly  proves  tliat  an  improvi-d  tesseral  harmonic  modi-l  of  the  earth 
IS  not  tile  way  to  go  in  order  to  improve  navigator  accuracy  due  to  the 
high  degree  of  the  model  requi ri-d  lo  give  even  modest  improvement. 

Ill  order  to  i-ompari-  ri-sulls  of  other  models  to  this  modi-l,  the-  correl- 
ation coi- 1 f ii- i eii t s tor  I’  and  g , and  for  F’  - 1’  and  g were 

I",  rr  f;v;  r 

computed  on  the  basis  of  the  tesseral  harmonic  model  with  random  phases, 
i'licse  correlation  coefficients  are  given  below  near  the  earth's  surface 
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li.  POINT  MASS  ANll  l.INK  MASS  MttDKl.S 


The  point  mass  and  line  mass  modi'Is  approximate  local  gravity 
anomalies  much  belter  than  the  tesseral  harmonic  model  which  has  primary 
use  in  describing  global  features. 

A simplified  analysis  of  a point  disturbance  begins  with  the  assiimvv- 
t ion  of  a point  mass  with  mass  m^| , 2-1)  x,  y 

coordinati'  frame  as  is  shown  in  l-'ig.  2. 


KIO.  2 IMINT  MASS  MODKl, 
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Vho  (1  ist  mbiuK'f  potfulial  tUu'  Id  at  a pi>uil  (x,  y)  is  simply 

V - p V p x"  » >■"  . Kxprt'ssiDus  for  tlu*  I'orct  s aiui  ^;ra(iitMUs  an* 

«»aslly  oblaliu’il  (>y  takluvJ  l-bo  rosportivo  partial  der  i vat  i vos , 
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whort*  r ^ x"  ♦ , 

riu’  rorrrl.'.t  ion  botw^uni  tin*  tour  (plant  ilios  in  (Kl)  can  bt*  found 
by  assuming  oaoli  of  t Iumii  i an  i ndopiMidcnl  int'asuroinont  , and  that  t lu* 
mi'asuromt’nl  s arc  obtained  alonp  .i  pri*scrlbcd  patli  in  I lie  x,  y plane, 
Ki>r  simplieily,  l lu*  path  ehosen  heiu'  is  I lu*  patli  v const  ant,  and  x 
proet't'ds  from  ^ at  a constant  rate  ( t lie  path  an  airplane  mif;lu 

follow),  flu*  informal  ii>n  matrix  is  then  ^:ivt'n  by 


L J 


and  t lu*  covarianct*  matrix  P is  simply  I , Normalized  off— diagonal 
elenu'uts  of  P r.lvt*  t lu*  corielatiou  cot*  f f i c i ent  s , A typical  inlc'{;ra— 
t ion  is  shown  below 
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^•.ravity  fail  i onu'f  t'v  .‘>ys(«'in  imp!  I’nu'iit  al  totis  . In  most  oas«'S,  1 lu'  f'.tavilN 
}’,  I'ail  i oiiu' ( rr  is  ustsi  as  a smsor  to  pit>vi^lo  tom'  |u'!*turlsit  ions  ilii*'  1 
j’.ravttv  anomaltos  that  t lu'  ai’oo  1 1' is»m«'t  oannot  pr»>vitU'  aiul  that  lontain 
t»it»  nuu’h  hl>;li  I' ro^pioiu*  n iiitormat  Ion  Ihir  an  a priori  rarth  t’.ra\Mlv  motitl 
to  appi'i'X  i ma  ( o . Instoait  of  prootssf  i nji  alon^:  f ho  path  iif  i n f <*p.  ra  f t Jif:  t h<' 
I’.railiont  inlormatlon  to  I’.lvo  tin*  l\>ro<'s  (ulili'h  loahs  t vi  proh|»'nr-  it  t lit' 
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with  t lu'  K'vt’s  a way  of  p locood  i tit*  directly  from  the  gradients 


tt>  the  forci's  without  augmenting  the  states  of  the  system.  In  fact, 
from  (Iti)  we  obtain 
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wlK're  y is  a correlation  liciglit  or  a mean  hciglit  of  the  vehicle  above 
the  disturbances.  A simple  gain  adjustment  is  all  that  is  neeiled  then 
to  provide  the  forci*  perturbations  from  t lu'  gradient  information. 
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I, INK  M;\SS  MtSiKI.S 


rite  line  mass  model  is  atiofher  possible  altertiati’  lo  locally  model- 
ing gravitational  perturbations.  I’o  di'termine  if  line  mass  nuidels  or 
point  mass  models  ri'sult  in  better  performance  will  require  survi'ving 
d.ita  to  discover  if  local  gravity  perturbations  are  more  accurately 
moiieli’d  with  line  masses  or  point  massi’s.  Surprisingly,  the  use  of 
point  masses  vs  1 I lu’  masses  jiroduces  rt'latively  minor  differences.  fh*’ 

1 1 lU'  mass  model  is  shown  in  Kig.  d. 
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The  disturbance  potential  V is  simply 
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pd  y, 

' 22  2^ - 
(x  +y  t7.  )‘ 


(18) 


where  now  ii  is  a mass  per  unit  li'tipth.  (Note  that  V in  (18)  Is 
independent  of  /..  ) Kxpressions  for  the  force  and  ^;rndient  expression  art' 
apain  obtained 
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I'lu'  cor  rt' 1 a I i on  between  the  various  elements  is  obtained  by  the 

procedure  described  in  (11)  and  ( 1 3 ) . The  Identical  substitutions  are 

made  here  as  in  the  point  mass  model,  the  I nt ep:rat i ons  carried  out,  and 

the  i-csults  are  that  Is  correlated  with  (and  with  nothing  else) 

with  correlation  coel' f ic  i ent  1 .'2  tl.  71  and  that  n - " is 

XX  ' vy 

correlated  with  1'^  (and  with  nothing  else)  with  correlation  coefl'icii'nl 
- l V -l).7l.  Although  the  correlation  Is  not  as  strong  as  with  t he 
point  mass  models,  it  Is  still  quite  high.  I’erhaps  it  is  appropriate 
in  the  final  model,  to  actually  use  a correlation  coefficient  between 
the  values  given  by  line  mass  models  and  point  mass  models.  Again, 
scale  heights  can  be  derived  to  relate  directly  the  gradient  perturba- 
tions with  the  force  perturbations. 


An  obvious  udvantanc'  of  those  lorrclntlon  imnlels  over  dynamtc  models 
of  the  earth's  gravity  gradient  perturbations,  deflections  of  the  vertical, 
and  perturbation  force  components,  Is  the  simplicity.  A dynamic  model 
involves  earth  surveying  to  "fit"  parameters  in  the  model,  which  may 
themselves  vary  substantially  over  different  regions  of  the  earth. 

An  added  advantage  is  tliat  the  correlation  models  do  not  involve  "aug- 
menting the  state"  of  the  system,  as  do  the  dynamic  models. 


Chapter  III 
BIAS  ESTIMATION 


Regardless  of  the  eventual  use  of  the  gravity  grad iometer , some  type 
of  signal  processing  will  be  required  to  minimize  the  effects  of  one  of 
the  largest  gravity  gradiometer  errors:  an  unknown  instrument  bias. 

As  will  be  discussed  in  Chapter  V,  a gradiometer  bias  has  the  same 
effect  as  a gyroscope  drift  when  used  in  an  inertial  navigator,  i.e., 
they  both  produce  unbounded  position  errors  with  increasing  time.  When 
the  gradiometer  is  used  as  the  sensor  to  perform  geodesy  experiments 
from  orbit,  the  gradiometer  bias  produces  a bias  in  the  estimate  of  the 
dominant  term  of  the  earth's  gravity  field. 

This  Chapter  contains  a study  of  several  attempts  to  eliminate  the 
gravity  gradiometer  bias  error. 

One  greatly  simplified  model  of  the  bias  estimation  problem  is 
given  below.  A single  (rotating)  gradiometer  supplies  two  pieces  of 
information  of  the  gravity  gradient  at  a point,  the  difference  of  two 
principal  elements  of  the  gravity  gradient  tensor  and  one  diagonal  coit>- 
ponent.  For  example. 
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22 
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12  • 


(20) 


A set  of  three  such  gradiometers  pointing  in  orthogonal  directions 
along  with  Laplace's  constraint  equation  (that  the  sum  of  the  diagonal 
elements  equal  zero)  is  sufficient  to  determine  the  gravity  gradient  at 
a point;  i.e.,  the  problem  below; 
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Page  bijuac 


is  invertible,  and  hence  solvable  for  x in  terms  of  z.  (It  is  not 

T - 1 T 

invertible  without  the  constraint  equation  Zy ) In  fact,  x = (H  H)  H z. 

Even  if  random  noise  v,  with  covariance  r.  is  inserted  into  each  of 

1 1 

the  first  six  equations  in  (21),  it  is  possible  to  obtain  a least - 

T ^1  —1  T —1 

squares  estimate  for  x in  terms  of  z.  In  fact,  x = (H  R H)  HR  z 

T -1  -1 

and  the  variance  of  the  estimate  error  P = (H  R H)  , Since  z^  is 
an  "exact"  measurement  though,  either  infinitesimal  values  for  the 
corresponding  element  of  R must  be  used  and  a limit  process  applies 
or  else  the  matrix  inversion  lemma  [10]  may  be  used.  The  matrix  inver- 
sion lemma  can  yield  the  improved  covariance  matrix  P from  the  covar- 
iance matrix  P' . P'  is  the  easily  computed  covariance  matrix  which 
does  not  include  the  beneficial  effects  of  the  exact  constraint  equation. 
The  formula  relating  P to  P'  is 

P = P'H^(H3”H^)”^HP' 

where  H reflects  the  constraint  equation  structure  that 
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The  advantage  of  the  matrix  inversion  lemma  procedure  is  that  it 
completely  avoids  the  limitlnit  processes  mentioned  earlier. 

Now  consider  the  problem  of  nu^;mentlng  the  state  vector  with  con- 
stant but  unknown  biases,  1 bins  per  iometer . Now 


’‘ii  - ’'22  " N 
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+ b„ 


’‘dd  - ’‘11  ^ •’d 
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Id 


4 b. 


^l  ' ’‘22  " ’‘dd  = « 


(2d) 


and  althou>;h  the  sum  of  ttie  biases  ’*  *^2  ^ ^d  obtainable  (from 

t impossible  to  obtain  any  other  combination  of 

the  biases  (in  particular,  the  individual  biases)  and,  tlicrefore,  it  is 
impossible  to  solve  for  x. 


The  only  useful  results  from  (2d)  are  obtained  if  there  is  some 
a priori  estimate  of  the  biases.  In  this  case,  the  usual  Kalman  filter 
techniques  will  allow  for  estimates  of  x to  be  obtained  from  noisy 
measurements,  however  the  conditioning  of  the  problem  is  of  course  de- 
pendent on  the  Initial  accuracy  of  b.  If  is  the  initial  informa- 
tion on  the  augmented  state  (P^^  ^0^^’  update  equations  are 

simply 


/s  ^ _ 1 T - 1 

’‘1  ^ ’‘0  ^ S ' 

1 , I + H'SrSl' 

1 (1 


(2d) 


where  II' 


is  the  measurement  matrix  for  the  augmented  state  vector. 


ami  K tlu'  covarianct'  ol  the  measurement  accuracy. 


riu'  aiivantane  of  llte  initial  information,  1 is  evident  in  the 

T -1  ‘ 

second  iijuat ion  ot  (2d).  Whereas  H'  K H'  is  not  invertible,  the  sum 

+ ir*H  *11'  IS  invertible  and  hence  the  Inverses  appearing  in  the  estim- 
ate update  equation  are  valid. 

More  sophisticated  schemes  for  estimating  the  biases  make  use  of 
dynamic  models  of  the  system,  us  opposed  to  the  model  in  which  both  the 
gradient  components  and  the  biases  were  constants.  Of  course,  since  an 
a prion  model  of  either  tlic  measurement  or  the  state  presumes  additional 
ini ormal ion  is  available,  it  should  be  expected  that  better  results  could 
be  obtaineii  at  tlie  expense  of  additional  complexity. 


Two  additional  bias  estimation  schemes  will  be  discussed.  Tlie  first 
makes  use  of  the  fact  that  the  bias  is  fixed  in  the  instrument  frame,  while 
what  the  gradiometer  measures  will  generally  be  fixed  in  an  inertial  frame. 

Ilotating  tlie  spin  axis  of  the  gradiometer  should  then  distinguish  between 
the  gradient  components  and  the  bias  terms.  The  second  method  again 
distinguishes  betwi'en  the  gradient  components  and  the  bias  terms  but  now, 
the  mechanism  is  via  a model  of  the  gradient  perturbations.  This  is  where 
much  of  the  mati'rial  from  Ch.  I can  be  used.  'M 

Heginning  first  with  the  rotation  scheme,  an  expression  for  the 

gradient  tensor  is  I 
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Assuming  the  gravity  gradiometer  spin  axis  is  rotated  about  the  1 axis 
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amt  H IS  lime  varying.  I >'  tact 


(2-n 
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A(;aiii  wi‘  have  the  cons  t ra  i at  equation,  + '^22  * "'13 

that  liir-  t O,  the  I 1 mi-- va r y 1 n^  H in  (24)  attri'es  with  the  H in 

(23). 

rile  staiularil  lest  for  observability  can  be  perforiiietl  on  the  matrix  II 
in  (24)  to  lietermine  the  maximal  rank  of  the  matrix 


It  turns  out  that  the  maximal  rank  of  is  7,  so  one  mode  is  still 

not  observable.  .Additional  al(;cbra  yields  the  result  that  the  mode 
b^  ~ 1 observable,  but  that  2 '3  b^  + ^ is  observable.  What 

this  says  is  that  it  is  still  not  possible  to  distinguish  between  and 

b.  This  result  is  easily  explainable  since  the  ^radiometer  spin  axis 
never  has  a vertical  component  with  the  sini;le  axis  rotation  scheme  we 
have  considered  here.  However,  if  a two  axis  rotation  scheme  where  the 
spin  axis  of  the  (^radiometer  spans  all  direct  ioii'^,  all  the  parameters  are 
observable,  and  hence  the  (gradient  components  and  individual  biases  ore 
obtainable  from  /. 

A problem  with  the  continuous  rotation  schemes  yet  to  be  discussed 
is  the  introduction  of  a (;radient  field  due  to  the  kinematics  of  the  r.  ' a- 
t ion  itself.  Tlie  vector  formula  for  this  induced  gradient  field  is 

d ^ 2 

r = [u'  \(..'  X r)]  uv  - I . (23) 


In  two  dimensions  the  t;radlent  field  reduces  to  simply 


If  the  spin 


axis  rotates  as  slowly  as  2t  50  sec,  the  required  accuracy  of  this 

rotation  rate  to  obtain  full  potential  of  the  fixed  instrument  accuraev 
-9 

is  s I part  in  10  Since  this  accuracy  is  not  feasible,  the  contin- 

uous rotation  schemes  are  ruled  out. 
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still,  some  use  can  be  made  of  the  preccdinK  analysis.  I'he  fact  that 
all  the  slates  atui  biases  ate  observable  with  tlie  continuous  two  axes 
rotation  sclieine  means  that  with  a sufficient  number  of  discrete  measure- 
ment points,  tile  same  information  will  be  available.  This  can  be  accom- 
plished in  several  ways.  One  scheme  is  to  "calibrate"  the  gravity  gradio- 
meter  at  a fixed  position,  in  several  orientations  to  obtain  the  bias.  II 
the  bias  is  truly  constant,  the  unit  could  then  be  used  at  new  locations 
with  the  predetermined  value  of  the  bias.  A second  scheme  is  to  simply 
repeat  measurements  with  a moving  base  gravity  gradiometer  with  the  instru- 
ment oriented  in  different  directions.  For  example,  in  a surveying  mission, 
either  with  an  airplane  or  satellite,  it  is  possible  to  retrace  a ground- 
track  while  making  measurements  in  (three)  different  orientations.  This 
determines  the  bias  and  stales  on-line,  in  which  case  the  effects  of  bias 
drift  coulil  be  minimized. 

The  second  scheme  for  distinguishing  the  gradient  components  from  the 
biases  is  via  a dynamic  model  of  the  gradient  componinit s . A simplifica- 
tion of  a model  to  appear  in  a later  chapter  appears  below.  Again,  letting 
X denote  the  gradient  components,  a linear  dynamic  model  of  x would  be 
X Fx  1 ’■'w  where  w is  a random  process.  Since  b 0,  there  is  again 
a mechanism  for  distinguishing  x and  b.  Adjoining  b to  x 
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llioiiph  II  K'll  ill  is  mil  i iivorl  i lilo  lor  iiso  In  I lio 

i'i|iiiil  loll,  I 111'  sum  of  I III'  I I'l  ms  on  I lio  rli^hl -liiinil  slilo  ol 
o vi'iil  111!  1 I y , II  will  111'  possllilo  lo  il  I s 1 i nun  1 sli  x I'rom 
usi'  ol'  ii  moili' I o!  I ho  oarlli  uriivlly  I'lolil  willi  iipplioii- 
nii  V I p III  I on  will  iippi'iir  Itilor, 
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I'liaptiM*  IV 


liKODKSY 


SrvtMiil  r N|U’ r i uuMit  s hav<‘  rmri>U'il  whlrli  liavt*  t lu’  tapabilily  oi  vt'iy 
ai'i  lira  t r I \ mrasurin};  I lu'  i;  ra  v i 1 a I i i>na  I lirltl  o{'  tin*  t'arth.  riu*s<*  };<'»>di'sy 
r\p<*r  i nu'Ml  s wliu’li  inrasuin*  t lu*  lupju*r  liarmonli's  of  t ho  t'aiMh  iiuluilt*  (1) 
lU}*,h“loNv  aiul  s a t r I I i t o- t iv- sa  1 1*  I M t t*  (S-S)  traoUtnj;;  (1!)  vt>inuor- 

or  1)  1 1 l M}*.  (T , O)  il  r a}*.—  1 ta***  sa  t o I I i t os  ; f d ) a 1 1 I mo  tor  m<*asuri*uu'M  t s f ri>m 
<>rhil;  f'l)  i>rhitii^;  p.ravity  i*.  r ad  t i>im‘ 1 1*  r nu'a.‘Oiia'mont  s ; (a)  t>t  lu'i'S  . 

riio  purpo^a*  i>f  till.**  oliaplor  i*.  to  ih*  1 1*  rm  1 no  t ho  aoouiaiy  1 1>  whuh  I In 
hipjioi*  harmonu's  of  t lu*  oartli'?'  ^;ravlty  fiold  ooiihl  ho  doiiiuu'd  with  an 
i>rh»tinf.;  f»ravity  rad  i omo  t I'l* , and  1 1>  I’ompari*  t lu'so  rosults  with  somi*  of 
t lu*  prov  ii>vi.s  I y montionod  «*\po r i inonl  . 

Tho  start  III}’  point  is  to  i'xpand  I lio  oartli's  f.iavity  fiolii  in  a a'lit’ 
of  tossoral  harmi>nli*s.  Koi'all  from  fl>  t ho  usual  roprosontat  ion  lU' 
t lu*  po I t urba t i on  poti'iitlal  as 
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It  Is  possiblo  to  dotluoo  analytii'  I'xprosstons  for  all  I lu*  }'.radit'nt  iomp<v- 

nonts,  but  for  simpliotty  aiul  sinoi*  it  will  lati'r  bo  assum«*d  that  t»nl\ 

r IS  moasuiaul,  an  oxprossu»n  lor  \'  will  suffuo  hoia*  (Aoluall\ 

rr  rr 

tlu*ro  aro  six  piooos  i>f  informal  li>n  availablo.) 
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All  Importaiil  point  to  l)i‘  noted  in  (28)  is  tiiat  although  tlie  out  ire 

Kfiivity  uradienl  tensor  at  a [loint  eannot  be  ree  oust  rue  ted  wltli  Just 

r , it  Is  possible  to  deduce  values  for  all  the  tesseral  harmonic  ci>- 

elficients  with  Just  r.^.  (or  in  (act,  with  any  single  tensor  component). 

In  fact,  the  measurement  i.  ’’  (r,  m.  \)  -t  v Is  linear  in  tlie  hurraonlc 

rr  ^ 

coefficients  (witli  wliite  noise  v added)  and  so  tlie  whole  theory  of 
linear  least-squares  estimation  theory  can  be  applied.  Ix.‘tting  the  state 
vector  X contain  all  tesseral  liarmonlc  coefficients, 

/.  llx  I V (29) 
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Kor  satellite  coverage  of  sufficiently  large  portions  of  the  earth,  the 
InfornKi t ion  matrix  (inverse  of  the  covariance  matrix)  can  bo  approximated 
as 

H^(r,tp,\^)  Hj(r,y,\,^) 
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(90) 


since  11  IS  a scalar;  (90)  makes  use  of  the  fact  that  the  satellite  is  in 
a polar  orbit,  and  hence,  latitude  coverage  is  continuous,  while  longitude 
c o ve rage  is  d I si‘ re t e . 
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Unless  111  111',  then  with  sufl'ieient  eovernt^e  1 ^ 0,  so  harmonic^  «llli 

il  i f t'erent  ni  decouple  and  the  only  lemaininn  problem  is 


(-2n  /1‘ 

\ (i-i)i«i-.Oti''iH(' ^2)1  — 


1.1 


whoro  T is  the  correlation  time  of  tlie  im*asurt'ment  noise  in  units  iif 
radians  of  orbit.  Rewriting;  (d2),  wt'  have 
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Unless  ( the  integral  in  (dd)  is  zero.  When  ('  (',  t lie  value 
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so  1 Is  diagonal  and  t lu'  I'leiiieiits  are 
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where  the  averam-  value  ol  sin  and  eos  Is  taken  to  be  *. 
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informat toil  matrix  for  the  C,  and  S has  diagonal  elements 

( m (in 
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This  gives  the  ^variance  ' for  a particular  or  of 
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Its  W/  n^(f+l)(e 


Again,  using  Kaula's  rule  of  thumb, 


^ X 10-*^ 


I 2it  RT 
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/r  \ 1 10® 

r hr  " TITI 


n = 3000  E 


T 10  sec  = 0.0116  rad 


R 0.1  E 


# orbit  = 3 months  ~ 1200  orbits. 


E’igure  '1  shows  a comparison  of  various  three-montli  geodesy  missions, 
including  the  current  knowledge,  counter  orbiting  satellites  (C.O,), 
gravity  gradiometers  (G.G.)  at  different  altitudes,  and  sa t el  1 i t e- 1 o- 
satellite  (S-f)  tracking. 
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HARMONIC  ORDER 


FIG.  -l  A CGMPAHISON  OF  VAHIOUS  TlUlKK-MOOTll  MISSIONS. 


Due  to  the  fact  that  the  gradiometer  yields  a derivative  measurement  j 

it  tends  to  amplify  the  high  frequeney  components  of  tlie  earth*s  gravity  j 

field.  This  is  evident  in  Fig.  4,  'I’he  gradiometer  outperforms  the  var- 
ious other  experiments  in  the  high  harmonic  range. 
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Chapter  V 


kAuian  filtering  for  gravity  gradiometer  augmenting 

INERTIAL  NAVIGATION  SYSTEMS 


INTRODUCTION 

Inertial  navigation  is  based  on  the  simple  principle  that  position 
is  given  by  the  double  integration  of  acceleration.  Accelerometers  do 
not  measure  accelerations  themselves  but  measure  specific  forces  which 
are  defined  as  all  forces  acting  per  unit  mass  with  the  exception  of  the 
gravity  force.  Hence,  in  order  to  obtain  true  accelerations,  the  specific 
force  due  to  the  gravity  must  be  subtracted  from  the  outputs  of  the  accel- 
erometers, in  other  words,  the  gravity  acceleration  must  be  added: 

O = f + g (38) 


where 

Q = acceleration  of  the  vehicle 

f = specific  force 

g = gravity  acceleration  of  the  earth. 

Current  inertial  navigation  systems  use  a "reference  ellipsoid 
model”  to  compute  the  gravity  acceleration  of  the  earth,  given  the  vehicle 
position.  Although  the  reference  ellipsoid  model  can  well  approximate 
the  real  gravity  field  of  the  earth,  the  difference  is  becoming  a major 
error  source  of  inertial  navigation  systems  because  of  rapid  hardware 
technology  advances  With  accelerometers  and  gyroscopes  [ll]. 

The  difference  between  the  actual  gravity  and  the  reference  ellip- 
soid model  may  be  expressed  in  terms  of  gravity  anomaly  (magnitude) 
and  del  lections  of  the  vertical  (angular  deviation),  [12] . 
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From  measurements  taken  at  12.5  nm  intervals  across  the  35th  parallel 
in  the  United  States,  the  standard  de\iation  and  the  correlation  distance 
of  the  vertical  deflection  were  deter’nined  as  5.2  arcseconds  and  25.1  nmi, 
respectively  [Gelb,  A. 1974] . The  worldwide  vertical  deflection  ensemble 
is  considered  to  have  8 arcsecond  rms  and  20  n mi  correlation  distance  [13]. 

Compensation  for  the  errors  causid  by  the  gravity  deflection  and  anom- 
aly is  one  of  the  principal  applications  of  gravity  gradiometers  currently 
under  development.  Gravity  gradiometers  measure  gravity  gradients  which 
are  related  to  gravity  acceleration  by  the  following  relation: 


where 

r = gravity  gradient  tensor  of  the  earth 

V = velocity  of  the  vehicle. 

Given  the  gravity  gradients  with  the  velocity  and  the  initial  valius 
of  the  gravity,  we  can  integrate  (39)  to  obtain  the  gravity.  However,  as 
is  well  known,  bias  errors  of  the  gravity  gradiometers  may  produce  unbounded 
position  errors  as  time  increases.  This  fact  is  easily  seen  by  the  follow- 
ing simple  example  illustrated  in  Kig.  5.  Assuming  that  only  error  source 
is  the  bias  of  the  gradiometer.  the  linearized  error  proj^agation  equation 
for  a single  horizontal  channel  with  constant  speed  may  be  given  by 


> 

II 

( lOa) 

Av  = Ag^ 

(10b) 

Ag  = vA  I’  + r Av 

(10c) 

X XX  XX 


where 

Ax  \v  and  Ag^  = estimation  errors  of  position,  velocity  and 
gravity  disturbance,  respectively 
b = M'  = bias  error  of  the  gradiometer. 

XX  '' 
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Fur  the  spherical  earth,  the  value  of  1'  is  constant  and  equal  to 
- ^ = -l-lOO  E,  where  R is  th»'  radius  of  the  earth.  Substituting  - 
for  r , we  can  inteerate  (40)  to  yield 

XX 


, vb  A 

Ax  = — t + — cos(u  t + k') 

' C 


\v  c=  — + A sin(.  t + v'') 
- s 

U.' 

s 


\g  c Au'  cos(..  t + 0) 
X s s 


where  v ind  v'  are  constants  deteriuined  by 


— + A sin  = initial  velocity  error 


A.  ^ cos  = initial  (gravity  error. 


Equation  (11a)  indicates  that  the  position  error  due  to  the  bias  of  the 
pradiometer  becomes  unbounded  with  increasing  time. 

In  order  to  overcome  this  difficulty.  Heller  [13]  proposed  a method  of 
using  gravity  gradiometer  as  an  external  aid  combined  with  a gravity  deflec- 
tion model.  He  obtained  a number  of  numerical  results  for  a single  horizon- 
tal channel,  assuming  velocity  reference  errors,  accelerometer  errors,  and 
gr.idiomoter  errors. 

In  the  following  sections,  we  try  to  obtain  an  analytic.il  solution  for 
Heller's  mechanization,  considering  the  gradiometer  error  as  the  only  error 
source.  Then,  we  extend  his  mechanization  to  estimate  the  bias  error  of  the 
gradiometer  (the  bias  error  in  this  case  means  the  difference  of  the  me.ins 
of  the  outputs  of  the  gradiometer  from  the  gravity  deflection  model). 


B.  SOLITION  FOR  A SINGLE  HORIZONTAI,  CHANNEL 


Til'  mech.inizat  ion  considered  here  is  the  same  as  " Oradiomet  er-as-,in- 
•xternal  aid"  (GAEA)  in  Heller  fid),  see  Fig.  (i.  The  gravity  obtained 
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FIG.  6 GRAinO'.KTKK  AS  AN  EXTKKNAL  NAVIGATION  AltJ  (GAEA)  [from  Heller  1975 


I 


by  the  rot'orenci-  ellipsoid  is  usid  I'or  luivigation  computation.  External 
velocity  information  is  provided  in  feedback  form  in  order  to  damp  the 
Schiller  oscillation.  Ciradiometer  measurements  are  combined  with  a gravity 
deflection  model  to  estimate  the  navigation  irrors  via  use  of  a Kalman  fil- 
ter. The  vehicle  spied  is  assumed  to  be  constant. 

Among  various  statistical  gravity  deflection  models  IHef.  1.)],  we 
choose  the  second  order  Markov  model  which  is  suitable  for  Kalman  filter 
implementation  because  the  governing  equations  may  be  written  in  tlie  fonii 
of  linear  differenti.il  equations  given  by 


tr?,  = -(-■'gf,  + Pgf, ' 
gP’  = 


(‘12) 


where 


f.  = vertical  deflection 
r, ' = aiigmi'iited  state 

P = Markov  p.irameter 

= zero  mean  wliite  noise  with  power  spectral  density  q^ 
and  q,  are  given  by 


= 


2,  146v 


q.  = 


•>  o 

•lK“f?ms 


(43) 


wh  re  1)  .\nd  ?,  are  the  correlation  distance  and  nns  of  the  vortical 

I'ms 

de  lection  respectively. 

In  this  section,  wo  consider  that  the  gradiometer  error  is  described 
by  a zero  mean  white  noise.  Tlien  the  measurement  equation  may  be  approxl- 
m.ited  by 


Z = v(r 

XX 


+ vV 


= -l^g^  + I'gf' 


+ vV 

g 


( l l.i) 
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whiTe 


i'  = output  of  Kradiomet  or 
XX 

(r  ) 1 computed  by  using  the  reference  ellipsoid  model 

XX  lit  XX 

= zero  mean  wliite  noist'  of  the  gradiometer  with  power  spectral 

density  V . 

K 


IS  given  by 


V = TAP 


T = averaging  time 


Ar  = rms  of  gradiometcr  error. 
K 


For  simplicity,  and  to  make  clear  the  effect  of  the  gravity  deflection,  we 
assume  that  tlic  only  position  error  source  is  due  to  gravity  deflection. 
Then,  tlie  error  propagation  equation  for  a single  horizontal  channel  may 
be  written  as 


Ax  = Av 


Av  = -a.'“  AX  - 2fu;  Av  + g. 


r = external  velocity  damping  coefficient 
A g 

0)  = ^ = Schuler  angular  frequency, 

s R 

As  is  clearly  seen,  ttie  gravity  deflection,  5,  and  its  augmented  state 
are  observable  by  the  gradiometer  measurement,  but  the  navigation  errors  in 
position  (Ax)  and  velocity  (Av)  are  not  observable.  Hence,  we  can  con- 
struct  Kalman  filter  for  f,  and  and  have  the  best  estimates  f,  and 

I'  given  by 


ge  = -&gf, + gi'  + K^(Z  + vPgP.  - vt3gp,') 


g?>=  -p3gf^  + Kp,  (Z  + vpgt  - vdgt’) 
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.ii't'  Kalin.in  filtrr 


wlu'ir  K aiui  K . 

f P * 

Now,  1 ho  osllination  orrt»r  in  gravity  (i  t s t ur  l)a]un‘ , «irivos 

t hr  orror  r(|uallon  (JTi)  , r.illior  i lian  tla*  t’ull  privity  d i s i nrhanoo , . 
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llu'  syninu'trtr  rmit  liu-iis  with  pii riiiiu' I < ' v .1  muy  lie  i i'.ulily  iliMwii  a-,  slusvii 
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wo  also  I'ind  I lio  sloady  slalo  Kalman  lillor  i;a  1 ns  (>ivon  by 
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K = 0 
r. 


(50) 


TIk’  cov  ir  i ;inci’  of  ( lu'  I’st  iiiiat  mii  itimi'  of  1 ho  <ioll''ction  ^ as 

(;iven  lay 
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F f 
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As  t ho  .icciiraoy  of 

oi-i'or  covai'ianco 

F = 8 ai'csoo, 
a'ms 


I ho 


I)  = 


Rraaiiomi't  oi'  impj'ovos  (M'  >0  oa' 

K 

lioca’oasos  moiaol  oai  i cal  ly  (Kit;.  ^ • 

DO  nma.  ami  v = 100  knots,  wo  hnv.' 


a > •■")  , 
Foa'  oxampli 
f com  (45) 


I ho 


foa- 


[i  = 10.75  (hr*') 

5 -3 

q,  =2.01  (aa.  rai  lia-  ) , 


For  f = 10  si-c,  .M'l;  = ' k’  ( wo  havo  fa'om  ( l ib) 

V = 4.07  \ 10”^(ha-~'S  . 


Hoiaco,  f t'oiii  ('ISO)  aiad  (51),  wo  havi' 


a = 15'1 

p"  = .Oils, 

F F 

I'hi'  ims  of  1 hi'  I'stimalion  oa'a'oa'  in  tho  vortical  iloflootaon  ilocroasos  to 
about  o lo  loiatti  of  tliat  witlioiit  tho  crail  i onio  I oa' . Ilowovi'a',  as  a > .• , 
oiai  claa  lacl  I'f  i s I i c root  appa'oactii's  tlio  oripiia  amt  tlio  otlii'r  appa'oaolaos 
inlinity.  Tli  a s fact  iaiilicatos  tlial  tlio  trfavity  itofloclioii  is  olitainoit  In 
a n I I'Ki'a  I i on  of  ttio  naoasiiriniont  , i|;iioa'int;  tlio  tlfnvity  ilofloction  nioito  1 . 
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ODAL  DESTABILIZATION 
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EIG.  H NOf'.V^LIZED  COVARIANCE  OF  VERTICAL  DEFIJ.CTION  ESTI'.UTION  ERROR. 


Subsl  itut  mu  (lie  ostim:itioii  orror  of  tlio  (gravity  deflect  ion  f-f  foe 

r in  (ISb)  and  condiictinn  tedious  calculation,  we  have  expressions  for 

the  cov-iriance  of  errors  in  position  (P  ) and  velocity  (P  ) given  by 
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covarianci's  of  I'l-rors  in  position  ami  velocity,  without 
grad i ometer  measurement 
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Among  numerical  examples  shown  in  KigJ^.  9,  t hi’  case  with  tlie  vehicle  veloc- 
ity loot)  knots  is  particularly  interesting  because  when  the  root  mean  S(|uari 
values  of  the  gravity  gradiomefer  error  are  around  10  K,  the  covariances  of 
position  and  velocity  errors  becimie  largi'r  than  tliose  witlioul  gradiome  t ers . 
This  fai’t  suggests  studying  the  power  spectral  density  of  the  estimation 
error  of  the  gravity  di’flection  which  is  given  by 


(PSD)  ~ 
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I'  4 ‘2*2 
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(53) 
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100  knots 


whtTf 


" \ U' 

Examining  (53)  and  Kig.  10a  wo  find  that  noar  7.oro  froqiuncy  tho  power  sih'C- 
t ral  density  of  the  istimation  error  of  the  gravity  deflection  is  always 
larger  (up  to  a factor  of  2)  than  This  means  that  even  if  tho  Kalman 

filter  gives  a smaller  covariance  as  shown  in  Kig.  8,  it  does  not  give  im- 
proved information  near  zero  frequency,  but  worse.  Hence,  when  the  power 
spi’Ctral  density  of  the  gravity  di'flection  error  at  Schuler  frequency  is 
l.irger  than  q.  , tin.'  covar i .inci’S  of  position  and  velocity  errors  are  larg- 
er with  gradiomet  I'rs  than  without. 

C.  BIAS  ESTINUTE 

In  this  section,  we  make  an  .ittempt  to  ext  I'lid  lU'ller's  mechanization 
to  estim.iti'  the  bias  error  of  the  gradiometer,  introducing  the  bias  as  an 
augminted  stati'  tn  the  Kalman  filter  discussed  in  the  previous  secti<in. 

Till'  bias  I'rror  in  this  case  means  the  difference  between  the  means  of  the 
outputs  of  tlie  gradiometers  and  the  gravity  deflection  model. 

The  measurement  equation  (■1-1)  may  be  rewrit  ti-n  as 

Z = - t^g^  + Pg^’-t  vb  + vVg  (5‘1) 

where  b is  the  bias  error  of  the  gradiometer  and  the  additional  state 
equation  is  simply 

b = 0 (55) 

The  full  system  consists  of  equations  (42),  (54)  and  (55).  Direct  appli- 

cation of  the  Kalman  filter  theory  to  this  problem  fails,  howiver.  Since 
the  bias  is  an  und i s t ur bable  mrl  neutrally  stable  mode,  the  Kalman  filtei 
gain  associati'd  with  this  mode  becomes  zero  after  I lu'  initial  uncertaintly 
disappears.  This  is  a typical  exampli-  of  Kalm.iti  filler  divergeiici.  Many 
cures  have  been  proposid  for  this  difficulty  such  .is  rest.irling,  minimum 
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FIG.  10a  POiVKR  .SPFXTRAI.  DENSITY  OF  EST 


100  knots 


poraii  SPECTRAL  DENSITY  OF  ESTIMATION  ERROR  OF  VERTICAL  DEFLECTION 
(Mf)flal  Destabilization). 


v.iri.iiice  obsiTVirs  with  i’ igeiiv-i  1 uo  constraints,  added  noise,  pole-shifting 
and  it'  stabilization  which  are  discussed  in  Bryson  [H].  Here,  we  use 
the  medal  di’s t abi  1 i za t ion  method  wliich  is  based  on  the  fact  that  the  steady- 
state  Kalman  filter  for  a system  with  an  (undisturbed)  unstable  mode  is 
stable.  As  the  amount  of  destabilization  increases,  the  absolute  value  of 
the  eigenvalue  associated  with  the  undisturbed  model  increases  from  zero. 
However,  thi-  covariance  of  tlie  estimation  errors  increases  too.  Some  num- 
erical results  .ire  showai  in  Fig.  8 .md  Table  1.  The  power  spectral  density 
of  the  estimation  I'rror  of  the  gravity  deflection  was  computed  numerically 
and  shown  in  Fig.  lOb.  Although  there  is  still  a hump  higher  than  the  orig- 
itial  value  , the  |ii>wer  spectr.il  density  of 
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Ch;ipter  VI 

KSTIMATION  OF  DISTRIBUTED  MASS  DENSITY 

Wi'  luivt'  been  trying  ii  completely  different  way  of  processing  gradio- 
meter  measurement  from  that  discussed  in  the  previous  chapter.  Instead 
t>f  equation  (39),  we  use  t lie  gradiometer  measurement  to  estimate  the  mass 
densit\-  distribution  of  the  earth,  then  compute  the  gravity  from  the  den- 
sity distribution.  This  approach  not  only  brings  up  a very  interesting 
problem,  i.e.,  filteiing  of  a distributed  system, but  also  seems  to  have  the 
following  practical  advantages: 

(a)  The  bias  of  the  gradiometer  does  not  cause  an  unbounded  position 
erixir  but,  at  most,  e.xcites  the  Schuler  oscillation.  This  is 
because  both  the  gravity  and  the  gravity  gradients  are  obtained 
by  spatial  integration  of  the  mass  density  and  time-integration 
of  the  gradiometer  measurement  is  not  needed. 

(b)  Sinci'  the  statistical  model  of  the  mass  density  distribution  is 
needed  only  at  the  boundary,  the  effect  of  the  model  error  is 
small.  Furthermore,  there  is  little  difficulty  in  extending  this 
•ipproach  to  the  actual  inhomogeneous  earth  density  field. 

On  the  other  hand,  obviously  this  method  requires  large  computei' 
capacity.  However,  since  tlie  correlation  distance  of  the  gravity  deflec- 
tion is  about  twenty  nautical  miles,  we  do  not  hav’e  to  estimate  the  den- 
sity distribution  over  large  areas  for  inertial  navigation  purposes.  This 
fact  may  relax  the  requirement  of  the  computer  capacity. 

So  far,  we  have  derived  the  partial  differential  equation  for  tlie 
density  distribution  and  obt.iined  the  filtering  algorithm,  including  dis- 
tributed Kalman  filter  g.iin.  Numerical  calculation  is  in  progress. 
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I’ONl'l.H.SlONS  AND  lil-X^OMMKNDAT IONS  I'OK  Dl'irrilKH  .STl'lA' 


Si'Vi'iMl  piMVity  nioili’lK,  i iicl  iidi  nt; 

1)  t('ssi‘r:»l  ha I'moiii c nunlcls 

2)  poinl  mass  aiul  line  mass  models 

;0  second  oi'der  random  pi'ocisss  models 

have  bi'en  st\Hlied.  As  a result  o!  I lie  fact  that  ini’rtial  navij;atoj-s  ari' 
pa  r t i on  1 ai' 1 \ si-nsilivi’  to  what  happi'iis  locally,  it  could  he  proved  ripoi- 
oiisly  that  Mi'  ti'ssei'al  harmonic  model  probably  could  lU'Vi-r  be  im|)lemenled 
line  to  the  larpie  number  of  parameters  lU'eded  to  obtain  a suitably  accui'ate 
local  p:ravtly  description. 

The  second  order  random  priu'i'ss  model  was  stiidieil  in  di'tail.  We  ob- 
tained an  analytical  solution  for  Heller's  iiu'chani'/.a  t i on  (.HAKA)  , nsinp  a 
more  simplifii'd  problem  formulation  which  still  retains  t hi’  part  I’ssential 
for  praii  i ome  t e r study.  I'he  solution  shows  that  the  covarianci'S  of  t lu- 
errors  in  position  and  velocity  are  larjiei-  than  those  without  (;i'ad  i ome  t er 
will'll  Schulei'  freiiuency  falls  within  the  bandwidth  of  the  estimation  error 
ill  the  vertical  deflection. 

We  extended  lleller'.s  Kalman  filter  to  estiniati'  the  bias  error  of  the 
rirail  i OHIO  t e r , usinr,  the  modal  des  I ab  i 1 i a 1 i on  mi'thod  to  avoid  Kalman  filter 
d 1 ve  rpi  nci' . 

I’oint  mass  and  line  mass  models  still  should  be  cons  i de  ri'd , due  to 
the  simplicity  with  which  they  could  be  i mpl  emeu  t I'd . I'Titure  woi'k  will 
contain  niiiiierical  resii'ts  based  on  these  model.;. 

Keciiitly,  .inotlii'i'  method  ot  obl.iininp  gravity  pei't  iirbat  i ons  has  coiiu' 
niidir  cons  i lie  r.i  t i on , t hi'  estimation  of  distributed  mass  di'ii.sity.  At  this 
point,  the  necessary  aii.ilytlcal  work  h.is  Ih'i'Ii  completed,  and  future  work 
will  include  numerical  results. 

Tlio  problem  of  bias  I'st  im.it  ion  ol  the  p:r.ivity  jir.id  i onie  t I'r  h.is  now 
been  stinlii'd  in  some  det.iil.  Afti'r  sevi'l'al  sclienies  weri'  Cons  i dei'i'd . the 


I 


most  f.ivoriiblo  motlmds  to  emorgo  are 


1)  If  possiblo,  retrace  a given  groundtrack  with  tlie  gradiometer 
in  tliree  different  orientations. 

2)  Calibrate  the  gravity  gradiometer  bias  at  a fixed  location  in 
many  different  orientations  to  obtain  the  bias  accurately  before 
the  instrument  is  used  as  a system  component. 

Finally,  the  results  of  using  a gravity  gradiometer  to  perform  a 
geodesy  mission  accurately  are  compared  with  competing  schemes.  A low 
orbiting  gradiometer  appears  to  bo  the  most  effective  way  of  obtaining 
tlK' tesseral  harmonics  of  the  earth  for  order  40  and  above,  i.e.,  the  high 
frequency  components.  In  the  range  below  the  40th  liarmonic,  the  geodesy 
mission  employing  counter-orbiting  drag  free  satellites  is  superior. 
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